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O . Abstract 

o. 

The tetra-maximal neutrino mixing pattern predicts a relatively large reactor mixing angle ^13 ~ 
8.4°, which is in good agreement with the latest best-fit value ^13 = 9°. However, its prediction 
r^ ' for 9i2 ~ 30.4° is inconsistent with current oscillation data at the 3cr C.L. We show that explicit 

o,. perturbations to the tetra-maximal mixing can naturally enhance 9i2 to its best-fit value 612 = 34°. 

^ '. 

r~| , Furthermore, we demonstrate that if the tetra-maximal mixing is produced by a certain flavor 



symmetry at a high-energy scale A = 10^^ GeV, signiflcant radiative corrections in the minimal 
super symmetric standard model can modify 612 to be compatible with experimental data at the 



(N 
> 

f~^ ! electroweak scale Aew = 10^ GeV. The predictions for ^13 ~ 8.4° and 623 = 45°, as well as the 



CP-violating phases p = a = —5 = 90°, are rather stable against radiative corrections. 
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I. INTRODUCTION 



Recent solar, atmospheric, reactor and accelerator neutrino experiments have provided us 
with compelling evidence that neutrinos are massive and lepton flavors are mixed [l|. The 
lepton flavor mixing is described by a 3 x 3 unitary matrix V, the Maki-Nakagawa-Sakata 
(MNS) matrix |2|, which can be parametrized by three rotation angles and three CP- violating 
phases. In the standard parametrization advocated by the Particle Data Group |3| and in 
Ref. y, the MNS matrix reads 

^ C12C13 S12C13 Srse-A fe'P 0^ 

V = -S12C23 - Ci2S23Si3e*'' C12C23 - Si2S23Sl3e*^ 
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S23C13 

S12S23 - Ci2C23Si3e'" -C12S23 - Sl2C23Sl3e'" C23C13 

where Sij = sin 6ij and c^- = cos 6ij (for ij = 12, 23, 13). The latest global analysis of current 
neutrino oscillation data yields 31° < 6112 < 37°, 36° < 6*23 < 53° and 4° < 6*13 < 13° at the 
3cr C.L., and the best-flt values of three mixing angles 612 = 34°, ^23 = 40° and ^13 = 9° |5|. 
Note that the global analysis has shown more than 3a evidence for a non-vanishing reactor 
mixing angle ^13 7^ 0, while the maximal atmospheric mixing ^23 = 45° is still allowed 
at the la C.L. However, three CP-violating phases {6,p,a) are entirely unconstrained. The 
smallest mixing angle ^13 and the Dirac CP-violating phase 6 will be measured in the ongoing 
and forthcoming neutrino oscillation experiments, while the Majorana CP-violating phases 
p and a can be constrained in the neutrinoless double beta decay experiments and colliders. 
How to understand the lepton flavor mixing pattern remains an open question in ele- 
mentary particle physics. Based on the observed neutrino mixing angles, however, several 
interesting constant mixing patterns have been proposed and widely discussed in the con- 
text of flavor symmetries. For instance, the tri-bimaximal mixing pattern with 612 ~ 35.3°, 
^23 = 45° and 6*13 = is in good agreement with current oscillation data [6|. Its predictions 
of ^23 = 45° and ^13 = have motivated a torrent of activities in the model building with dis- 
crete flavor symmetries, which give rise to the tri-bimaximal neutrino mixing at the leading 
order (seeje.g., [7| and references therein). Nevertheless, the latest result from the T2K ex- 
periment |8|, in which z/^ — )■ v^. oscillations have been observed, indicates 5.0° < 6*13 < 16.0° 
for the normal mass hierarchy and 5.8° < ^13 < 17.8° for the inverted mass hierarchy at the 
90% C.L. The best-flt value from T2K data is 6*13 ^ 10°, which is also consistent with the 
global-flt analysis. Moreover, the MINOS experiment has recently reported the observation 



of Ufj, — !■ Ue oscillations, disfavoring the assumption of 6^13 = at the 1.5a C.L. Possible 
ways to realize a relatively large ^13 have recently been discussed in Ref. 9|. Considering 
the experimental tendency of a non- vanishing or even relatively-large ^13, we argue that the 
tetra-maximal mixing pattern 10 1 
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with 9i2 ~ 30.4°, 623 = 45°, and 613 ~ 8.4° may serve as a better starting point to search 
for the true symmetry underlying the lepton flavor mixing. The prediction of ^13 ~ 8.4° 
from the tetra-maximal mixing pattern is in excellent agreement with the latest neutrino 
oscillation data, while that of 623 = 45° and the maximal CP- violating phase 6 = —90° may 
hint at a certain flavor symmetry. 

In the present work, we consider possible deviations from the tetra-maximal mixing pat- 
tern. The motivation for such an investigation is two- fold. First, the solar mixing angle 
predicted by the tetra-maximal mixing pattern is 612 ~ 30.4°, which is much smaller than 
the best-fit value 6*12 = 34° and even not lying in the 3a range from the global analysis. Sec- 
ond, if the tetra-maximal mixing pattern is obtained by assuming a certain flavor symmetry, 
which in general works at a high-energy scale, the mixing angles will receive significant ra- 
diative corrections when running from the symmetry scale to the low-energy scale, at which 
the mixing angles are actually measured in oscillation experiments. One immediate ques- 
tion is whether it is possible to increase 6*12 to the observed value by introducing explicit 
perturbations to the tetra-maximal mixing pattern or by taking into account the radiative 
corrections, while both 6'23 and 613 remain consistent with experimental data. We have found 
the answer is affirmative. 

The remaining part of this work is organized as follows. In Sec. II, we show that only one 
perturbation parameter is enough to increase 612 to its best-fit value, and leads to an even 
larger ^13. Furthermore, we demonstrate that both the maximal atmospheric mixing angle 
6'23 = 45° and the maximal CP-violating phase 6 = —90° are not affected by the perturbation 
to any order. In Sec. Ill, we explicitly solve the renormalization group equations (RGB's) for 
neutrino mixing parameters within the minimal supersymmetric standard model (MSSM), 



and confirm that the observed 612 can be obtained at the electroweak scale Aew = 100 GeV 
from the tetra-maximal mixing pattern at the high-energy scale A = lO^'^ GeV. On the 
other hand, both ^13 and 623, as well as three CP-violating phases (5, p, a), are rather stable 
against the radiative corrections. Finally, we summarize our conclusions in Sec. IV. 

II. EXPLICIT PERTURBATIONS 

First of all, we briefly recall the tetra-maximal mixing pattern and its salient features. 



As shown in Ref. 10|, the tetra-maximal mixing matrix in Eq. (2) can be decomposed into 



four maximal rotations 

V = Pi® 023(vr/4, 7r/2) ® Oi3(vr/4, 0) ® Oi2(7r/4, 0) ® Oi3(7r/4, vr) , (3) 

where Pi = Diag{l, l,i}, and Oij{6ij,6ij) is a rotation with the angle 6ij and the phase 6ij 
in the complex i-j plane for ij = 12,23, 13. Comparing between Eq. (1) and Eq. (2), one 
can extract three neutrino mixing angles 

tan^i2 = 2-y2, tan^as = 1 , sin^i3 = -(2 - ^2) , (4) 

or explicitly 612 ~ 30.4°, 6'23 = 45°, 6*13 ^ 8.4°. Furthermore, three CP-violating phases are 
p = a = —6 = 90° can be obtained by redefining the phases of charged-lepton fields and 
recasting the tetra-maximal mixing pattern into the standard form in Eq. (1). 

Although 6*13 ^ 8.4° is relatively large and well compatible with the recent T2K and 
MINOS results, 612 ~ 30.4° is much smaller than the best-fit value 612 ~ 34°, and even not 
covered in the 3a range, i.e., 31° < 612 < 37°. In order to increase 612 significantly but 
not change much both 6'23 and 613 such that all three mixing angles become consistent with 
current oscillation data, we have to slightly modify the tetra-maximal mixing pattern in a 
proper manner. One straightforward way is just to introduce a small perturbation to the 
maximal rotation in the 1-2 complex plane. In this case, the MNS matrix turns out to be 

V = Pi® 023(7r/4, 7r/2) ® 013(71/ A, 0) ® Oi2(7r/4 + £12, 0) ® Oi3(7r/4, n) , (5) 

with £12 ^ 1. In order to see how the neutrino mixing angles are changed, we expand the 



MNS matrix in Eq. (5) with respect to the perturbation parameter to the second order 

f -V2 2 v^ ^ 
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Some discussions are in order: 



1. Comparing Eq. (6) with the standard parametrization in Eq. (1), one can immediately 
derive three neutrino mixing angles 
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and tan 6*23 = 1. Taking the best-fit value of 612 = 34°, we find £12 ~ 0.11 from the 
first identity in Eq. (7). Inserting £12 ~ 0.11 into the second identity, one can obtain 
6*13 ^ 10.8°, which is in good agreement with the T2K result 
upper limit 6*13 < 13° from the global- fit analysis |5|. 



and close to the 3a 



2. Note that the maximal atmospheric mixing angle, i.e., 6'23 = 45°, is not modified, 
because there is a generalized permutation symmetry in the tetra-maximal mixing 
matrix in Eq. (3) and its perturbed version in Eq. (5): ^23^^ = U* for both U = V 
and U = V, where P23 denotes the exchange of the second and third rows. Such a 
permutation symmetry originates from the pure phase matrix Pi and the first maximal 
rotation 023(7r/4, 7r/2), which are the only complex matrices. Hence we have tan6'23 = 
l^sl/lK-sl = 1 that is independent of the perturbation. 

3. The tetra-maximal mixing pattern predicts the maximal CP-violating phase 6 = —90° 
and the Jarlskog invariant JT" = Im [V'e2V"^3V'g3V?2] = —1/32 for leptonic CP viola- 



tion 



10| . Taking account of the perturbation, we arrive at 



^ 32 
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Therefore, the magnitude of leptonic CP violation, measured by \J\, is enhanced due 
to a positive £12, which has been implemented to increase both 612 and ^13. However, 
we can prove that the maximal CP-violating phase 5 = —90° is maintained to any 
order of perturbations. The proof is as follows: 

(a) Since the matrix elements Vei (for i = 1,2,3) are always real, the Jarlskog in- 
variant is i7 = Im [V'e2V'^3Vg3V'^*2] = Ve2K3lHi [^/r2^M3] • ^u the other hand, we 
have J' = si2Ci2S23C23Si3cf3 sin 5 in the standard parametrization, or equivalently 
i7 = |V'ei||K;2||K3| sin(5/2, where tan6'23 = 1 has been input. Therefore, the 
CP-violating phase is determined by 

2v;2v;3im [v^^v^s] 
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(b) The unitarity of the MNS matrix V leads to the normalization condition jV^ip + 



IT/ |2 I IT/ |2 
I yfi2 1 + I V3 1 



1 and the orthogonality condition V^^^ + V^2 + ^/^3 



latter condition is guaranteed by P23V^ = V*. Thus we obtain 
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Subtracting Eq. (10) from Eq. (11) and using another normalization relation 
I Kip -I- 2|V^ip = 1, one can verify \Vei\ = 2 |lm [V^*2V^3] |, which together with 
Eq. (9) leads to | sin5| = 1. As long as the perturbations are small, the sign of J' 
is determined by the tetra-maximal mixing matrix V and thus 6 = —90° is valid 
to any order of perturbations. 

4. In the basis where the charged-lepton mass matrix is diagonal, we can reconstruct the 
neutrino mass matrix by the MNS matrix V and neutrino masses rrii (for i = 1, 2, 3) 
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With the help of the identity P23K = y*, it is straightforward to show that 
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where P23 = -^23^ = -^23 and m = Diagjmi, m2, ma}. As a consequence of Eq. (13), 
we can get 

where M^^ denotes the matrix element of M^ (for a, /3 = e, /i, r). Such a special struc- 
ture of My, which can give rise to both maximal atmospheric mixing 6*23 = 45° and 
maximal CP- violating phase 5 = ±90°, may result from a certain flavor symmetry. 
Since both V and the tetra-maximal mixing matrix V share the same permutation 
symmetry, the perturbation under consideration does not spoil the symmetry of neu- 
trino mass matrix in Eq. (14). 

If we generalize the perturbation scheme in Eq. (5) and add small perturbations to all 
four maximal rotations in Eq. (3), the MNS matrix can be written as 
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where the higher-order terms 0(£^) have been neglected. In this case, the neutrino mix- 
ing angles receive corrections from all these perturbations. Comparing Eq. (15) with the 
standard parametrization in Eq. (1), we obtain 
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to the first order of perturbations. Since £13 contributes only to the solar mixing angle 612, we 
can switch off both £12 and £23, and choose £13 ~ —0.11 to obtain the best-fit value of 6*12 = 



34°. The unique and efficient way to enhance both 612 and 613 is to adjust the perturbation 
parameter £12 as we have done before. Now that £23 breaks the permutation symmetry 
-^23^ = V^*, it induces deviation of ^23 from the maximal mixing, and also invalidates the 
maximal CP-violating phase. In addition, the symmetry relations for the reconstructed 
neutrino mass matrix in Eq. (14) are not respected. 

III. RADIATIVE CORRECTIONS 

Now we proceed to consider another possible deviation from the tetra-maximal mixing 
pattern, i.e., the renormalization group (RG) corrections. As already mentioned in Sec. I, 
the flavor symmetries generating the tetra-maximal mixing are generally preserved at high- 
energy scales, such as the grand unification scale (e.g., Agut = 10^^ GeV) or the hypothetical 
seesaw scale (e.g., A = 10^^ GeV), while the neutrino mixing parameters are determined or 
constrained in neutrino oscillation experiments at low-energy scales. The gap between the 
high-energy predictions and the low-energy measurements is bridged by the RG evolution. 
The RG running effects may then serve as an explanation for the discrepancy between the 
flavor symmetric mixing pattern and observables. 

The RGB's for neutrino mixing parameters have been derived within various theoretical 
frameworks ll|. In the supersymmetric theories with a large tan/3, it has been found 
that the RG evolution may lead to significant modifications to the mixing parameters, in 
particular the solar mixing angle Ou (see e.g., Ref. [12] and references therein). To be 
explicit, we write down the RGB's for neutrino mixing angles in the MSSM in the leading- 
order approximation |l3l ]. 
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where 9ij = d^/dt with t = \n{fi/fio), ( = Am'^^J Aml^ with Am^Qj = ml — mf ^ 
7.6 X 10"^ eV^ and |Am^tj^| = |m| — mg] ~ 2.3 x 10"^ eV^ |5[, and i/r denotes the Yukawa 
coupling of tau charged lepton. Note that the terms of 0(6*13) have been safely neglected in 
Bq. (17). In addition, we have defined C2{p-~a) = cos2(p — a), C2p = cos2p and so on. Some 

8 



general features of the RGE's are summarized as follows: 

• Since \Am?^^^\/ ^m1^^ ^ 30, one can observe that 612 in general receives more remark- 
able RG corrections than 623 and ^13. Furthermore, when running from a high-energy 
scale A = 10^^ GeV to the electroweak scale Aew = 100 GeV, the radiative correc- 
tions to 6^12 are always positive, i.e., 6^12 (A) < 6^12 (Aew)) which is independent of the 
neutrino mass hierarchy and the CP-violating phases. Therefore, if the tri-bimaximal 
mixing pattern with 612 ~ 35.3° is assumed at the high-energy scale, the RG effects 
will lead to an even larger 612 at the low-energy scale, which may be in conflict with 
current oscillation data [ij]. Note that 6*12 receives negative corrections in the stan- 
dard model, but there are no visible RG effects in this case due to the absence of tan /3 
enhancement. 

• As for 6^23 and 6^13, the RG corrections could be either positive or negative, depending 
on the neutrino mass hierarchies. More explicitly, we have 6'23 < in the normal 
hierarchy (NH) case with m^ > ?Ti2 > frii, while 623 > in the inverted hierarchy (IH) 
case with m2 > rrii > m^. The sign of ^13 further depends on the CP-violating phases. 

On the other hand, the RGE's for the CP-violating phases are approximately given by [l3| 

Vr i (2 2 n"^1'5i2-52p+ (l + C)"^2C?2S2a , mx'm2c\2s\j^S2{p-„) 
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SttH ^^^^ '" AmL^(l + C) ' AmL 

where the terms of 0(6*13) have been ignored. For the tetra- maximal mixing pattern, one 
can insert the initial condition p = a = —6 = 90° into Eq. (18), and immediately obtain 
p = (T = at the leading order, indicating that all the Majorana phases are rather stable 
against RG corrections. As for 5, we have 5 oc -, — "^3(m2-mi) — ^ -^ ^^^ ^ neutrino mass 
hierarchies. Therefore, there is no enhancement factor boosting the RG running of 5, and 
(5 = —90° is also stabilized at any energy scales. This has also been confirmed in our 
numerical calculations. 

To illustrate the RG corrections to the tetra-maximal mixing pattern, we have numerically 
solved the full set of RGE's for neutrino mixing angles. Some comments are in order: 
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FIG. 1: The RG evolution of Qx2 in the MSSM with tan/3 = 15. The sohd and dashed Unes 
correspond to r7ii(A) = and mi (A) = 0.05 eV in the NH case, or 1113 (A) = and 771,3 (A) = 0.05 eV 
in the IH case, where A = 10^^ GeV is a typical seesaw scale. The colored bands indicate the RG 
corrections to the mixing angle if the lightest neutrino mass is varying in the range of (0 ^ 0.05) eV. 
The allowed ranges of Q\-2 from the global analysis [5| are also shown as shaded areas. 

1. In Fig. 1, the RG evolution of ^12 in the MSSM with tan/3 = 15 is shown. Note that 
we have assumed the tetr a- maximal mixing pattern at a cutoff scale A = 10^^ GeV 
(i.e., potentially the seesaw scale), and allowed the lightest neutrino mass ?tii in the 
NH case (or 7713 in the IH case) to vary in the range (0 ~ 0.05) eV. As seen from 
Fig. 1, although the initial value of Qyi at the cutoff scale deviates more than 3cr from 
its best-fit value, the RG effects can enhance Q\2 in a very efficient way so as to fit 
the experimental data. In the NH case, Q\2 can be perfectly consistent with the low- 
scale measurements for the chosen tan/3 and neutrino masses. In the IH case, the RG 
corrections with tan/3 = 15 seem to be too large. This behavior can be understood 
by noting that the first identity in Eq. (17) reduces to Q\2 oc {m2 + m\)l(m2 — mi) in 
the limit of p = cr = 90°. In the NH case we have (7712 + mi)/{m2 — nii) > 1, whereas 
in the IH case (7722 + 77ii)/(7722 — fni) > 10^ is expected. Therefore, a small tan/3 (e.g., 
tan /3 < 20) is required in the IH case to avoid the overlarge RG corrections. 

2. In Fig. 2, we show the RG evolution of 6'23 and ^13 in the MSSM with tan/3 = 15. 
Compared with the case of 612, the evolution of both ^23 and ^13 is insignificant due to 
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FIG. 2: The RG evolution of ^13 and ^23 in the MSSM with tan/3 = 15. The sohd and dashed hnes 
correspond to mi (A) = and ?n,i(A) = 0.05 eV in the NH case, or ?n,3(A) = and 7713 (A) = 0.05 eV 
in the IH case, where A = 10^^ GeV is a typical seesaw scale. The colored bands indicate the RG 
corrections to the mixing angles if the lightest neutrino mass is varying in the range of (0 ~ 0.05) eV. 

the suppression from !\vp?^^yj\!\m}^^^. The radiative corrections to Q21, and Q\z cannot 
exceed 0.1° no matter whether the NH or IH is assumed. It is worth mentioning that 
^13 = at the leading order, which can be seen by inserting p = o = —b = 90° into 
Eq. (18). The mild evolution of 6*13 in Fig. 2 is actually attributed to the high-order 
terms in the RGB's. Hence the predictions of 623 = 45° and ^13 ^ 8.4° from the 
tetra-maximal mixing pattern are rather stable against radiative corrections. 

3. One may wonder if we can acquire significant RG corrections to 623 and 613 by assuming 
a sizable tan/3 (i.e., tan/3 > 20) or a nearly degenerate mass spectrum (i.e., m^ > 
0.2 eV). Unfortunately, this is impossible, because it will cause too large corrections 
to 612- In Fig. 3, we depict the allowed regions in the plane of tan/3 and the lightest 
neutrino mass mi in the NH case (or 7713 in the IH case) by requiring 6*12 (Aew) to be in 
its 1(7, 2cr and 3a ranges. As shown in Fig. 3, both tan/3 and the lightest neutrino mass 
are severely constrained. In the NH case, tan/3 can be relatively large if mi < 0.05 eV, 
whereas tan/3 < 20 in the IH case for the whole range of 7713. In the nearly degenerate 
limit m,i ^ m,2 ~ m,3, we need a rather small tan/3, which may lead to tensions with 
some supersymmetric models. 
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FIG. 3: Allowed regions of tan/3 and rrii in the NH case (ms in the IH case), where the tetra- 
maxinial mixing pattern is assumed at A = 10"*^^ GeV and the 9i2 at Aew = 10^ GeV is required 
to be in the la, 2a and 3cj ranges. The shaded regions refer to the NH case, while the regions 
between lines to the IH case. 

Finally, we stress that our discussions on the RGB's stick to the effective theory approach 
featuring the tendency of generality. In a class of seesaw models with flavor symmetries, the 
heavy seesaw particles may possess non-degenerate masses. In this case, the RG running 
between seesaw thresholds should be considered accordingly |l5l]. In general, as long as 
one works in supersymmetric models, the threshold effects should not be significant because 
of the non-renormalization theorem. In the non-supersymmetric models, there might be 
large threshold corrections to neutrino mixing angles due to the existence of the Higgs 
self-coupling. However, the details of the threshold behavior depend on the specific flavor 
structure of the model, and in principle should be considered properly in model buildings. 



IV. CONCLUSIONS 



Recently, the long-baseline accelerator experiments T2K and MINOS have observed u^ — )■ 
Uf, oscillations, indicating a relatively large 6'i3, which will be soon tested in Double CHOOZ 
and Daya Bay reactor neutrino experiments. The latest global analysis of current oscillation 
data points to a best-fit value of ^13 = 9° [5| . If such a relatively-large 613 is confirmed in the 
near future, our understanding of the neutrino mixing pattern through flavor symmetries 
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would be dramatically changed. So far, various flavor symmetries have been implemented 
to generate the tri-bimaximal mixing pattern, which predicts 612 ~ 35.3°, 623 = 45° and 
6'i3 = 0. Now that a vanishing ^13 is disfavored with a more than 3a significance |5[, we are 
well motivated to consider significant corrections to the tri-bimaximal mixing or to search 
for another constant mixing pattern with a relatively large ^13. 



In this paper, we concentrate on the so-called tetra-maximal mixing pattern [lOJ, which 
predicts 6*13 ^ 8.4° together with 6*12 ~ 30.4° and 623 = 45°. Although a relatively large 
^13 is predicted, the solar mixing angle 612 seems too small to be consistent with current 
oscillation data. We propose two feasible ways to solve this problem. First, since the tetra- 
maximal mixing pattern can be written as a product of four maximal rotations, one may 
introduce explicit perturbations to one or more maximal rotation angles. We demonstrate 
that only one perturbation parameter is enough to enhance 612 to its best-fit value, and the 
maximal atmospheric mixing angle 623 = 45° and maximal CP- violating phase 6 = —90° are 
maintained to any order of perturbations. Second, if the tetra-maximal mixing pattern is 
produced by a certain flavor symmetry at a high-energy scale A = 10^^ GeV, the radiative 
corrections governed by the RGB's can successfully enhance 612 to its best-fit value at the 
electroweak scale Aew = 10^ GeV. We explicitly show that this is really the case in the 
MSSM by solving the full set of RGB's. In addition, ^13 and 623, as well as three CP-violating 
phases, are found to be rather stable against the radiative corrections. 

It is worthwhile to remark that the flavor symmetry underlying the tetra-maximal mixing 
pattern deserves further studies. The predictions for neutrino mixing parameters from the 
tetra-maximal mixing pattern with or without corrections will be soon tested in a number 
of precision neutrino oscillation experiments. 
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